The Violation of Stokes-Einstein Relation in Supercooled water 
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By confining water in nanopores, so narrow that the liquid cannot freeze, it is possible to explore 
its properties well below its homogeneous nucleation temperature Th~ 235 K. In particular, the 
dynamical parameters of water can be measured down to 180 K approaching the suggested glass 
transition temperature T 9 « 165 K. Here we present experimental evidence, obtained from Nuclear 
Magnetic Resonance and Quasi-Elastic Neutron Scattering spectroscopies, of a well defined decou- 
pling of transport properties (the self-diffusion coefficient and the average translational relaxation 
time), which implies the breakdown of the Stokes-Einstein relation. We further show that such a 
non-monotonic decoupling reflects the characteristics of the recently observed dynamic crossover, 
at about 225 K, between the two dynamical behaviors known as fragile and strong, which is a 
consequence of a change in the hydrogen bond structure of liquid water. 

PACS numbers: 61.20.Lc, 61. 12. Ex, 64.70.Pf 



Despite its fundamental importance in science and 
technology, the physical properties of water are far from 
being completely understood. The liquid state of wa- 
ter is unusual especially at low temperatures 0, 0, 0- 
For example, contrary to other liquids, water behaves 
as if there exists a singular temperature toward which 
its thermodynamical properties, such as compressibility, 
thermal expansion coefficient, and specific heat, diverge 
PJ . The efforts of scientists from many disciplines to seek 
a coherent explanation for this unusual behavior, in com- 
bination with its wide range of impacts, make water one 
of the most important open questions in science today. 
On the other hand, the nature of the glass transition 
(GT) of water represents another challenging subject for 
current research 4]. Dynamical measurements in glass 
forming liquids have shown a dramatic slowdown of both 
macroscopic (viscosity T] and self-diffusion coefficient D) 
and microscopic (average translational correlation time 
t) observables, as temperature is lowered towards the 
GT temperature T g . Accordingly, a comprehension of 
the GT has been sought through the study of the dy- 
namics at the molecular level, which, despite all efforts, 
has not yet been completely understood Keep- 
ing in mind the "complexities" of both low-temperature 
water and its GT, we present here direct measurements 
of two dynamical parameters of water: the self-diffusion 
coefficient and the average translational relaxation time, 
in the temperature range from 280 to 190 K, obtained by 
Nuclear Magnetic Resonance (NMR) and Quasi-Elastic 
Neutron Scattering (QENS) experiments, respectively. 

Bulk water can be supercooled below its melting 
temperature (Tm) down to « 235 K, below which it 
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inevitably crystallizes; it is just in such supercooled 
metastable state that the anomalies in its thermodynam- 
ical properties are most pronounced, showing a power law 
divergence towards a singular temperature Tg = 228 K. 
At ambient pressure, water can exist in a glassy form be- 
low 135 K. Depending on T and P, glassy water has two 
amorphous phases with different structures: a low (LDA) 
and a high (HDA) density amorphous ice; thus it shows 
a polymorphism. LDA can be formed from HDA and 
vice versa; LDA, if heated, undergoes a glass-to-liquid 
transition transforming into a highly viscous fluid, then 
crystallizes into cubic ice at Tx ~ 150 K. Thus an ex- 
perimentally inaccessible T region exists in bulk water 
between Th and T\ ■ Experiments performed within this 
interval could be of fundamental interest for understand- 
ing the many open questions on the physics of water. 
For example, the presence of a first order liquid-liquid 
transition line (LLTL), the precise location of its T g , 
recently suggested at about 165 K and the ex- 

istence of a fragile-to-strong dynamic crossover (FSC) 
on approaching T g from the liquid side [T(ij ]. The exis- 
tence of a LLTL leads to conjecture that liquid water 
possesses a low-temperature second critical point (pre- 
dicted to be located at T c w 220 K,P C w 1 Kbar) @], 
below which it can switch from one phase, a high den- 
sity liquid (HDL), to another phase, a low density liquid 
(LDL), whose corresponding vitreous forms arc the HDA 
and LDA, respectively. The difference between the two 
liquid phases lies in the water structure: in the HDL, 
the local tetrahedrally coordinated hydrogen-bond net- 
work is not fully developed, whereas in the LDL a more 
open, lo cally ice-like, hydrogen-bond network is fully de- 
veloped Thus, near T c , water is a mixture of both 
LDL and HDL phases associated with a diverging den- 
sity fluctuation. At higher temperatures, the two liquid 
phases are indistinguishable. Lowering temperature or 
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increasing pressure will result in an increase of the LDL 
phase with respect to the HDL phase. The FSC can be 
identified by the temperature at which transport proper- 
ties, like the shear viscosity r\ or the inverse self-diffusion 
coefficient l/D, cross over from a non-Arrhenius (fragile) 
to an Arrhcnius (strong) behavior on approaching T g . 

A possibility to enter this inaccessible temperature 
range of water, named "no-man's-land" , is now shown 
by confining water in nano-size pores 0, 0, 0, Il5| . 
When contained within these pores, water does not crys- 
tallize, and can be supercooled well below Tjj- Vycor 
pores (a porous hydrophilic silica glass), micellar 

systems or layered vermiculite clay |T^ are examples of 
confining nano-structures. The latter systems have been 
used to explore the Arrhcnius behavior of the dielectric 
relaxation time (td) of very deeply supercooled water. 

The FSC was recently confirmed by a QENS experi- 
ment, which measured the T and P dependences of the 
average translational relaxation time (tt) for water con- 
fined in nanopores of silica glass |l6l Il7|. In particular, 
as the temperature is lowered, a (tt) versus 1/T plot 
exhibits a cusp-like crossover from a non-Arrhenius to 
an Arrhcnius behavior at a temperature Tj,(P). This 
crossover temperature decreases steadily upon increas- 
ing P, until it intersects the Tr line of bulk water at 
P ~ 1.6 Kbar. Beyond this point, the FSC can no longer 
be identified. These results, suggestive of the existence 
of the two liquid phases, have been explained in a molec- 
ular dynamics (MD) simulation study by considering the 
existence of a critical point. The MD study shows that 
the FSC line coincides with the line of specific heat max- 
ima C™ ax , called Widom line. The Widom line is the 
critical isochore above the critical point in the one-phase 
region [lsf . Moreover, it is observed that crossing this 
line corresponds to a change in the T dependence of the 
dynamics. More precisely, the calculated water diffusion 
coefficient, D(T), changes according to a FSC, while the 
structural and thermodynamic properties change from 
those of HDL to those of LDL. 



Results and Discussion 

In this report, we present a detailed study done using 
two different experimental techniques, neutron scatter- 
ing and NMR, to probe dynamical properties of confined 
water at low temperatures, well inside the inaccessible 
region of bulk water. Our main aim is to clarify the 
properties of water as a glass forming material, measur- 
ing directly with NMR spectroscopy, the self diffusion 
coefficient I? as a function of temperature, and compar- 
ing the obtained results with the translational relaxation 
time (tt) measured by QENS. (tt) is a quantity propor- 
tional to the viscosity ij. These new measurements enable 
one to compare the proportionality of the transport co- 
efficient l/D and viscosity rj, and can provide a test for 
the theoretical description of dynamics at the molecular 
scale of glass forming materials. According to Ito et al. 



[To| . the FSC can be intimately connected to the pres- 
ence of a thermodynamic event in liquid water: i.e. the 
temperature dependence of the inverse self-diffusion coef- 
ficient does not follow that of viscosity or inverse mutual 
diffusion coefficient. In other words, supercooled water 
must show, on approaching T g , the marked decoupling of 
translational diffusion coefficient from viscosity or rota- 
tional correlation time, as recently observed in some su- 
percooled liquids |2fl l2ll |22| . Here we confirm by means 
of NMR data the existence of a FSC in supercooled wa- 
ter, as proposed by the QENS and MD studies [Hll7LIT| | 
and show, for the first time, that the Stokes-Einstein re- 
lation (SER) breaks down in different ways on both the 
fragile and strong sides of FSC. 

Figure 1A shows a log-linear plot of the inverse of the 
self-diffusion coefficient of water 1 /D measured by NMR 
as a function of 1/T for the fully hydrated MCM-41-S 
samples with pore diameters of 14 A and 18 A. Figure 
IB reports the average translational relaxation time (tt), 
obtained by analyzing QENS spectra of the same sam- 
ples, versus 1/T, u sing the Relaxing-Cage Model (RCM) 
as shown in Ref. [l(| E3 • As it can be observed from 
Figure 1A and IB, the measured values of D and (tt) 
are independent of the pore size of the samples. This 
indicates that NMR field-gradient measurements, hav- 
ing a length scale larger than the size of the pores, are 
insensitive to the system geometry. In Fig. 1A, the 
solid line denotes the fit of the data to a Vogel-Fulcher- 
Tamman (VFT) law l/D = l/D n exp(BT /(T - T )), 
where 1/D a = 2.4- 10 7 (s/m 2 ), B = 1.775, and T = 187 
K. B is a constant providing a measure of the system 
fragility and To the ideal glass transition temperature. 
The short dotted line denotes the fit to an Arrhcnius 
law l/D = 1 / Dqexp(Ea / ksT) , where we keep the same 
1/D value as in the VFT law fit, and E A = 3.98 
Kcal/mol. Figure IB shows the (tt) data at ambient 
pressure. The dashed lines denote the VFT law fit, 
and the dotted lines the Arrhenius law fit, with the 
same pre-factor tq. We obtained the following values: 
E A = 5.4 Kcal/mol, T = 200 K. The consequence 
of insisting on the same pre-factor in both the VFT 
and the Arrhenius laws results in an equation determin- 
ing the crossover temperature Tl in the following form: 
1/T L = l/T -Bk B /E A . We obtained T L = 224.5 K from 
the l/D data and T L = 225.8 K from the (t t ) data. The 
agreement between NMR and QENS results is thus sat- 
isfactory, especially regarding the two relevant quantities 
E A and Tl- Figure 1A and IB also show that within the 
size range of 10 to 20 A, the crossover temperature Tl 
is independent of pore size. As previously mentioned, a 
FSC occurring at 228 K has been proposed by Ito et al. 

for water, which is fragile at room and moderately 
supercooled temperatures but near the glass transition 
temperature it is shown to be a strong liquid by dielec- 
tric relaxation measurements [T^ . The interpretation of 
this transition as a variant of the structural arrest transi- 
tion (as predicted by the ideal mode coupling theory) was 
the essence of the recent QENS study of the structural 
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FIG. 1: Figure 1A shows, for the fully hydrated MCM-41- 
S samples with diameter 14 A and 18 A, the inverse of the 
self-diffusion coefficient of water D measured by NMR as a 
function of 1/T in a log-linear scale. The solid line denotes 
the fit of the data to a Vogel-Fulcher-Tamman (VFT) relation. 
The short dotted line denotes the fit to an Arrhenius law with 
the same prefactor 1/Do- Figure IB reports, as a function of 
1/T, the average translational relaxation time (tt) obtained 
from QENS spectra in the same experimental conditions of 
the NMR experiment. The dashed line denotes the VFT law 
fit, and the dotted line the Arrhenius law fit with the same 
prefactor To- In both panels, the values of fitting parameters 
are shown. 



relaxation time and MD study of the self-diffusion coef- 
ficient EE 12 EH These NMR results presented above 
constitute, by means of a direct measurement of the self- 
diffusion coefficient of supercooled water, an independent 
confirmation of the existence of FSC in water. 

Let us now focus on the SER that relates the self- 
diffusion coefficient D, viscosity rj, and temperature T as 
D (x T/rj, which, as it is well known, is usually accurate 
for normal and high temperature liquids. Since (tt) is 
proportional to the viscosity, we examine the relationship 
between D and (tt) in figure 2. In this figure the quan- 
tity D (tt) /T is reported as a function of T. Dots and 
squares represent its values coming from the experimen- 
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FIG. 2: The figure reports the quantity D (tt) /T as a func- 
tion of T. Dots and squares represent its values coming from 
the experimental data of D and (tt) in samples with diam- 
eter 14 A and 18 A, respectively. The dotted line represents 
same quantity obtained by using the fitting values obtained 
from the data reported in Fig 1A and IB. 



tal data of samples with pore diameters of 14 and 18 A, 
respectively, whereas the dotted line represents the same 
quantity obtained using the corresponding fitted lines 
reported in figure 1A and IB. The temperature depen- 
dence of D (tt) /T shows that this quantity is constant at 
higher T, but increases steeply as T goes below the FSC 
temperature. Furthermore, it shows a small bump at the 
FSC temperature, in accordance with the predictions of 
a recent theoretical study E3- Therefore, in the super- 
cooled region the temperature behavior of D and (tt) 
is inconsistent with SER, signalling a marked decoupling 
between these two transport parameters, on decreasing 
T . In recent experimental studies on some supercooled 
liquids, it has been reported that SER breaks down as the 
glass transition is approached. The self-diffusion coeffi- 
cient shows, as far as water in the present experiment is 
concerned, an enhancement of orders of ma gnit ude from 
what expected from SER [U H |H I2I l2^|. These 
decouplings of the transport coefficients, observed as a 
SER violation, have been attributed to the occurrence 
of dynamical heterogeneities in structural glass formers 
[20t I22I I25I l26| . Thus, in supercooled liquids there exist 
regions of varying dynamics, i.e. fluctuations that domi- 
nate their transport properties near the glass transition. 
The extent of such decouplings may depend on the mate- 
rial and the microscopic details of the specific transport 
parameters. 

The SER breakdown can be described using a scal- 
ing concept, in particular, the law D ~ t - ^, where 
£ = a(T)/ f3(T) with a and /3 being temperature depen- 
dent scaling exponents of D and r, respectively [27|]. Re- 
cently, it has been shown that for tris-naphthylbenzene 
(a fragile glass former) £ = 0.77 [22] , whereas an MD 
simulation of Lennard-Jones binary mixture has given 
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FSC. This supports the hypothesis that liquid water is 
consisting of a mixture of two different liquid structures 
(the LDL and HDL phases) . Accordingly, a liquid- liquid 
phase separation line exists in the P — T plane with a 
liquid- liquid critical point as its end point. Remarkably, 
we give the first experimental proof of the existence of 
a violation of SER above and below the FSC in water, 
i.e. in both the fragile and strong regimes of supercooled 
water. This clearly reflects the decoupling of transport 
coefficients of the liquid when temperature is lowered to- 
ward T g . This latter result certainly constitutes a new 
element which serves to clarify one of the most intriguing 
properties of water. 



FIG. 3: The scaling plot in a log-log scale of D vs (tt)- Red 
dots are data corresponding to temperatures above Tl, i.e. 
when water is in the fragile glass phase, whereas blues dots 
correspond to the strong Arrhenius region. Two different scal- 
ing behaviors exist above and below the temperature of the 
FST. In the fragile region the scaling exponent is £ ~ 0.74 
(dotted line) and ~ 2/3 in the strong side (solid line). Dashed 
line represent the situation in which the SER holds, D ~ r _1 . 



£ = 0.75 |2^. By using such an approach, we will discuss 
our SER results for confined supercooled water. Figure 3 
shows the D vs (tt) plot in a log-log scale. The red dots 
represent data corresponding to temperatures above Tl, 
where water behaves as a fragile glass former, and the 
blue dots pertain to the strong Arrhenius region. As it 
can be observed, the data clearly show two different scal- 
ing behaviors above and below the FSC temperature, in 
particular £ ~ 0.74 on the fragile side (dotted line) and 
~ 2/3 on the strong side (solid line). Dashed line repre- 
sents the situation in which SER holds, D ~ r _1 . These 
results agree with those obtained in tris-naphthylbenzene 
[23 and, more specifically, with those of a recent theo- 
retical study in which the decoupling of transport coef- 
ficients in supercooled liquids was investigated by using 
two class of models, one describing diffusion in a strong 
glass former, and the other in a fragile one [2^ . The 
main result of this study is that, while in the fragile case 
the SER violation is weakly dependent on the dimension- 
ality d, with £ = 0.73, in the strong case the violation is 
sensitive to d, going as D ~ t~ 2 / 3 for d = 1, and as 
D ~ r~ - 95 for d = 3. On considering the geometry of 
the system that we have used in our experiment to con- 
fine water (ID cylindrical tubes, with a length of some 
fim and pore diameters of 14 A and 18 A), the scaling 
plot shown in Figure 3 compares remarkably well with 
the findings of the theoretical investigation [23 , on both 
fragile and strong sides. 

In summary, we explore dynamical properties of water 
in a deeply supercooled regime (well inside the "no-man's 
land") by means of NMR and QENS experiments, which 
separately give a conclusive proof of the existence of a 



Methods 

Water was confined in micellar templated mesoporous 
silica matrices MCM-41-S, which have 1-D cylindrical 
pores with a length of some fim arranged in 2-D hexag- 
onal arrays, synthesized following a similar method for 
the previous synthesis of MCM-48-S 29]. The MCM-41- 
S materials are the same as those previously used in the 
QENS study of confined water [13. Pore size was de- 
termined using nitrogen absorption-desorption technique 
[l6L Investigated samples have hydration levels of 

H ~ 0.5 (0.5 gram H 2 per gram of MCM-41-S), ob- 
tained by exposing dry powder samples to water vapor at 
room temperature in a closed chamber. This water con- 
fining system can be regarded as one of the most suitable 
adsorbent models currently available |30L l3l| . 

High-resolution QENS spectroscopy method was used 
to determine the temperature dependence of (tt) for 
confined water. Because neutrons can easily penetrate 
the wall of sample cell and because they are predomi- 
nantly scattered by hydrogen atoms in water, rather than 
by the matrices containing them, incoherent QENS is 
an appropriate tool for our study. Using two separate 
high-resolution QENS spectrometers, we were able to 
measure the translational-relaxation time from 0.2 ps to 
10,000 ps over the whole temperature range under study. 
The experiments were performed at both the High-Flux 
Backscattering (HFBS) and the Disc-Chopper Time-of- 
Flight (DCS) spectrometers in the NIST Center for Neu- 
tron Research (NIST NCNR). The two spectrometers 
have two widely different dynamic ranges (for the cho- 
sen experimental setup), one with an energy resolution 
of 0.8 AieV (HFBS) and a dynamic range of ± 11 fieV H3, 
and the other with an energy resolution of 20 fieV (DCS) 
and a dynamic range of ± 0.5 meV 33] in order to be 
able to extract the broad range of relaxation times from 
the measured spectra. The experiment was done at a sc- 
ries of temperatures, covering both below and above the 
transition temperature, and the data were analyzed using 
Relaxing-Cage Model to extract the average translational 
relaxation time (tt). 

The NMR experiments on fully hydrated MCM-41-S 
samples with pore diameters of 18 and 14 A, were per- 
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formed at ambient pressure using a Bruker AVANCE 
NMR spectrometer, operating at 700 MHz proton res- 
onance frequency. The self-diffusion coefficient of water 
D was measured with the pulsed gradient spin-echo tech- 
nique (PGSE) jMEi, in the temperature range 190 K 
- 298 K (with an accuracy of ± 0.2 K). The T depen- 
dence of the chemical shift of methanol was used as a 
T standard. All details about the NMR experiment and 
the sample properties are reported elsewhere 0. The 
reported D values were derived from the measured mean 
square displacement (r 2 (i)) of molecules diffusing along 
the NMR pulsed field gradients direction r, during the 
time interval t. 
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